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Muon Bremsstrahlung and Muonic Pair Production in Air Showers. 
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The objective of this work is to report on the modifications in air shower development due to muon 
bremsstrahlung and muonic pair production. In order to do that we have implemented new muon bremsstrahlung 
and muonic pair production procedures in the AIRES air shower simulation system, and have used it to simulate 
ultra high energy showers in different conditions. 

The influence of the mentioned processes in the global development of the air shower is important for primary 
particles of large zenith angles, while they do not introduce significant changes in the position of the shower 
maximum. 



1. Introduction 

The physics of ultra high energy cosmic rays 
plays an important role in our days. These 
primary particles can not be detected directly. 
When an ultra high energy astroparticle interacts 
with an atom of the Earth's atmosphere, it pro- 
duces a number of secondary particles that con- 
tinue interacting and generating more secondary 
particles. This process is generally known as the 
air shower. 

We have developed a set of programs to sim- 
ulate air showers and manage all the output 
data. Such simulating system is known as 
AIRES (AIRshower Extended Simulations) @] . 

We started working on the topic of the elec- 
tromagnetic processes in air showers some years 
ago. We analyzed the modification in the shower 
development due to the reduction of the electron 
bremsstrahlung and electron pair production by 
the LPM effect and the dielectric suppression [|| . 
We also studied the influence of the geomagnetic 
field in an air shower ||. 

The main goal of this work is to analyze 
two processes that take place in an air shower: 
muon bremsstrahlung and muonic pair produc- 
tion (electron and positron). At energies high 



"Fellow of CONICET (Argentina) 



enough, these processes become more important 
than ionization, and therefore these mechanisms 
account for virtually all the energy losses for high 
energy muons. 

In order to study the modifications that these 
effects introduce in an air shower, we developed 
new procedures for such mechanisms and incor- 
porated them in the AIRES air shower simulation 
system jl]. The AIRES code has then been used 
as a realistic air shower simulator to generate the 
data used to make our analysis. 

This work is organized as follows: we start in 
section 2 with a brief summary of the theory of 
the muon bremsstrahlung and muonic pair pro- 
duction; in section 3 we show the results of our 
simulations; finally we present our analysis and 
comments in the conclusion section. 



2. Theory 

2.1. Theory of the Muon Bremsstrahlung 

The cross section for muon bremsstrahlung 
(MB) is calculated by the standard method of 
QED Q| similary as in the case of electron 
bremsstrahlung . 

The first approach to the MB theory was due to 
Bcthe and Heitler |5|-|7j . They considered in their 
calculation the screening of the atomic electrons. 
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After this first formulation some corrections 
were introduced. Kelner, Kokoulin and Petrukhin 
H also included the bremsstrahlung with the 
atomic electrons. The nuclear form factor was 
investigated by Christy and Kusaka Q for the 
first time and then by Erlykin [fTof . Petrukhin 
and Shestakov jll]] found that the influence of 
the nuclear form factor is more important than 
the ones predicted by the previous papers. This 
last results have been confirmed by Andreev et 
al. Q who also considered the excitation of the 
nucleus. 

2.2. Theory of Muonic Pair Production 

In the lowest significant order of perturbation 
theory the muonic pair production(MPP) is a 
process of 4th order in QED. 

Racah jl3j was the first to calculate the MPP 
cross section in the relativistic region without 
taking into accuont the atomic and the nuclear 
form factor. Thereafter, Kelner Jl4| included 
the correction due to the screening of the atomic 
electrons. The analytical expression for any de- 
gree of screening was introduced by Kokoulin and 
Petrukhin fl5| . Those authors also took into ac- 
count (T^| the correction due to the nuclear form 
factor. We wish to emphasize that the influence 
of the nuclear size is more important when the en- 
ergy tranferred to the pair is large Jl6| . This last 
case is important for the air showers and there- 
fore the nuclear size effect needs to be taken into 
account in the simulations. 

3. Air Showers Simulations 

We have calculated the total cross section, and 
equivalently, the mean free path (MFP), for both 
MB and MPP processes. 

Figure [j] shows the MFP in g/cm 2 for MB, 
MPP and muon-nucleon interaction as functions 
of the kinetic energy of the initial muon. The in- 
formation in this figure can be compared with the 
depth of the atmosphere. The vertical depth of 
the atmosphere is about 1000 g/cm 2 while in the 
case of very inclined showers (85° of zenith angle) 
is about 9000 g/cm 2 . Therefore, the probabilities 
of MB and MPP will not be very large except 
in the case of large zenith angles. The muonic 




Figure 1. Mean free path: MB, MPP and muon 
nuclear interaction vs the kinetic energy of the 
initial muon. 



component of the showers (number of muons di- 
vided by the number of electrons and gammas) at 
ground level becomes very important for zenith 
angles larger than 60° (see, for example, figure 2 
in reference ||); this is another reason to expect 
that the influence of both effects will be more ap- 
preciable under those conditions. 

The MFP's for MB and MPP diminish when 
the kinetic energy of the muon increases and 
therefore, the influence of both process will be 
more important for large energies. 

Due to the fact that the MFP of muon-nucleus 
interaction ^ is larger (more than one order of 
magnitude) than the corresponding for MB we 
do not take this effect into account in our calcu- 
lations. 

In order to analyze the influence of MB and 
MPP in the development of air showers initi- 
ated by ultra high astroparticles we have per- 
formed simulations using the AIRES program Q 
with different initial conditions: primary particles 
(protons, irons, muons), primary energies (10 18 , 
10 19 and 3 x 10 20 eV) and zenith angles. 

2 We used the same parameterization of the total cross sec- 
tion used by GEANT|17| 
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To start with the analysis of the simulated 
data, let us consider the case of a single muon 
(eventually produced during the development of 
a given shower) . This particle may generate a sec- 
ondary shower if the processes of MB and MPP 
are taken into account. This effect is clearly illus- 
trated in figure ||a where the longitudinal devel- 
opment of all charged particles is plotted versus 
the slant depth. In this case the initial particle is 
a muon of I0 14 eV. When MB and MPP are not 
taken into account there is practically no produc- 
tion of new particles. On the other hand, when 
the effects are considered, it is possible to appre- 
ciate a secondary shower. 
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Figure 2. a) Longitudinal development of all 
charged particles vs the slant depth. Pri- 
mary particle (Energy): muon (10 14 eV). 
b)Longitudinal development of muon fractional 
energy vs the slant depth. 



In figure ||b the longitudinal development of the 
muon fractional energy (muon energy/primary 
energy) is plotted. When the effects are not taken 
into account the muon practically does not loss 
energy during all its path, while if the effects are 
considered, the muon energy loss is significant. 

We have also studied the modifications in the 
global observables of the shower, for example 



when the primary particle is a proton. 

We have plotted in figures || and || different 
obsevables for 3 x 10 20 eV proton shower with 
zenith angle 85°, taking and not taking into a 
account the effects of MB and MPP. Figure || 
represents the longitudinal development of elec- 
trons and positrons (e + / _ ) versus the slant depth. 
The number of e +// ~ increase when the effects are 
taken into account. The differences between the 
two cases are clearly noticeable in the tail of the 
showers. The longitudinal development in energy 




MB-MPP On 



MB-MPP Off 



1500 3000 4500 6000 7500 9000 



X (g/cm 2 ) 



Figure 3. Longitudinal development of electrons 
and positrons vs the slant dept. Primary particle 
(Energy) proton (3 x 10 20 eV). 



of muons is plotted in figure [|. In this case the 
difference is about 15% at the tail. The energy of 
muons diminishes when the effects are considered. 

For the longitudinal development of muons ver- 
sus the slant depth the difference between the 
cases taking and not taking into account the ef- 
fects is less significant, about 2 % (not plotted 
here). 

The modification in the shower development 
due to the MB and MPP diminishes for small 
zenith angle of the primary particle. We did not 
find significant modification of the shower devel- 
opment for zenith angles smaller than 45°. 
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4. Conclusions 

Due to MB and MPP, a high energy muon in 
the shower can generate a secondary shower. This 
practically does not occur if the mentioned mech- 
anisms are not taken into account. 
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We are performing more simulations in order 
to make a complete study of the influence of MB 
and MPP in an air shower [fL8[ . 
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